The Fulde-Ferrell-Larkin-Ovchinnikov states for the d-wave superconductor in the 

two-dimensional orthorhombic lattice 
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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state of a two-dimensional (2D) orthorhombic lat- 
tice superconductor is studied based on the Bogoliubov-de-Gennes equations. It is illustrated that 
the 2D FFLO state is suppressed and only one-dimensional (ID) stripe state is stable. The stripe 
changes its orientation with the increasing Zeeman field. There exists a crossover region where the 
gap structure has some local 2D features. These results are significantly different from those of 
the tetragonal lattice system. The local density of states is also studied which can be checked and 
compared with experiments in future. 
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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, 
known as the finite momentum of the Copper pair, was 
predicted in the mid of It will occur when the 

Pauli paramagnetism effect dominates over the orbital 
effect 3 . Thus the low dimensional layered superconduct- 
ing (SC) materials with the magnetic field parallel to the 
SC layers are the strong candidates to realize the FFLO 
state. Recently, signatures for possible FFLO state 
have been found in various crystal systems, such as the 
heavy fermion materials CeCoIns^— the organic super- 
conductors A - (BETS) 2 GaCl 4 ^ A - (BETS) 2 FeCl^, 
k - (BEDT - TTF) 2 Cu(NCS) 2 r^ and the iron-based 
materials LiFeAs^ On the other hand, the realiza- 
tion of FFLO state in the cold atom system is studied 
intensivelyi 3 - - — and recently the experimental signatures 
are found in the one-dimensional optical latticed 

The experimental development for the FFLO state in 
low dimensional systems has attracted much attention. 
Theoretically the calculation based on the lattice model 
is of great interest and the results may compare with 
the experimental results of the crystal systems or the 
cold atom system. One of the fundamental issues of the 
FFLO state is the detailed gap structure, which can be 
studied through minimizing the free energy^, or through 
the self-consistent calculation based on the Eilenberger 
cquationj 20 i 21 At the mean time, the Bogoliubov-de- 
Gennes (BdG) technique has been a powerful tool to 
study various imhomogenous SC states and it can also 
obtain the SC order parameter self-consistently. Thus it 
is also an effective tool to study the FFLO state espe- 
cially in lattice system. Previously, based on the BdG 
equations, the FFLO state was studied intensively in the 
tetragonal lattice system. Many groups have reported 
their results and numerically the gap structure depends 
on the pairing symmetry^ - — For s-wave pairing sym- 
metry, the one-dimensional (ID) stripe-like pattern was 
rcported i 22 ' 23 For d-wave paring symmetry, recently it 
was proposed that a transition from the ID stripe-like 
pattern to the two-dimensional (2D) checkerboard pat- 



tern will occur as the exchange field increases^ And the 
2D pattern may change to the ID pattern in presence of 
the impuritiesj 24 ! 25 

The above results of the gap structure for the tetrago- 
nal lattice system are significantly different from the pre- 
vious theoretical results in the isotropic systemS ) 27 i 28 as 
discussed by Ref. [26], indicating that the intrinsic sym- 
metry of the crystal lattice should play an important role 
in the gap structure of the FFLO states. For a real SC 
material, there often exists a structural transition from 
tetragonal- to orthorhombic-lattice with the variation of 
the doping level. Some possible microscopic orders, such 
as the stripe order— or the nematic order— ~— may also 
lead to a weak anisotropy in the ab-plane. Furthermore, 
the ID FFLO stripe in the tetragonal lattice itself may 
induce the in-plane anisotropy. Therefore, the studies of 
the FFLO states in an orthorhombic lattice system is of 
great interest while so far little attention is paid to the 
FFLO state on this system. 

In this paper, motivated by the above consideration, 
we calculate the spatially distributed order parameter 
self-consistently based on the BdG equations on a 2D 
lattice with x — y anisotropic hopping. Our results show 
that the 2D solution of the SC gap structure will be sup- 
pressed and disappears as the orthorhombicity strength 
7 = t x /t y increases. For the case of 7 = 1.05, only stable 
ID stripe state exists and it is also found that the stripe 
changes its orientation with the increasing Zeeman field. 
A crossover region exists between these two stripe states 
and some local 2D features are obtained in this region. 
The local density of states (LDOS) is also studied to dis- 
tinguish the different states. 

We start from the model on the 2D lattice with the 
Zeeman splitting effect. The model Hamiltonian can be 
written as, 
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where tij are the hopping constants and /i is the chemical 
potential, ah is the Zeeman energy term, caused by the 
in-plane magnetic field, with a = ±1 representing the 
spin-up and spin-down electrons respectively. 

This Hamiltonian can be diagonalized by solving the 
BdG equations, 
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where is expressed by 
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The SC order parameter and the local electron density 
rii are obtained self-consistently: 
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Here /(x) is the Fermi distribution function. Vy is the 
pairing strength. In the present work, we consider the 
nearest neighbor (NN) pairing with Vij — V5 i: j±a. For 
the tetragonal lattice system the NN pairing will re- 
produce the <i a: 2_ y 2-pairing symmetry. The d-wave or- 
der parameter at the site i can be defined as Af = 



1/4(A, 



i-\-x 



A, 



Ai,i+y ~ Ai,i-y). For the case of 



orthorhombic lattice the four-fold symmetry is broken so 
that a s-wave component will be induced and is expected 
to increase as the orthorhombicity strength 7 increases. 
The s-wave component is defined as A? = l/4(A ii+x + 



A, 



Ai_i- y ). Since we consider only weak 



x — y isotropy in the present work with 7 < 1.05. Our nu- 
merical results show that the s-wave component is quite 
small [A 3 / A d < 0.03]. Thus in our following presented 
results we neglect the small s-wave component and use 
the above d-wave order parameter as the definition of the 
on-site order parameter. We also define the The magne- 
tization mi as rrii — (Sf) — h/2{n^ — n^}. 
The LDOS is expressed by 



[\u^6(E n -u) + \v%\ 2 6(E n +u>)}, (6) 



where the delta function 8{x) has been approximated by 

t/tt(x 2 + r 2 ) with r = 0.02. 

In the following calculation, we consider the NN hop- 
ping with the hopping constant in x direction being 
t x = 1.0 and that in y direction t v = t x /j. Here, 7 
represents the orthorhombicity strength. The pairing po- 
tential V and the filling electron density n are chosen as 
V = 1.3 and n = 0.84 (hole-doped samples with doping 
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FIG. 1: Phase diagram of an 2D orthorhombic lattice su- 
perconductor in the parallel magnetic field, h = q^bH is 
the exchange field and 7 the orthorhombicity strength. The 
crossover region is indicated by the shadow. 



6 = 0.16). The calculation is made on a 48 x 48 lattice 
with the periodic boundary condition, and the randomly 
distributed initial values of the order parameters are cho- 
sen. The 10 x 10 supercell is used to calculate the LDOS. 

Our main results are summarized in Fig. 1, where 
the phase diagram is plotted. As one can see, for low 
anisotropy, the whole SC state is divided into three re- 
gions, namely, the uniform SC state, the ID FFLO state, 
and the 2D FFLO state, which is consistent with the 
previous results^ For 7 = 1.05 and at low fields, the 
SC state still includes the uniform SC state and the ID 
FFLO state. At 0.22 < h < 0.23, as indicated by the 
shadow, there exists a crossover region. In this region, 
the gap structure forms the coexistence of the ID stripe- 
like pattern and some local 2D features. As the exchange 
field increases further, as seen, another ID FFLO state 
shows up. The orientation of the stripe pattern is dif- 
ferent from that of the low field, namely, for the lower 
field FFLO state the stripe is parallel to the x direction 
(stripe I) and that of the higher field one is parallel to 
the y direction (stripe II). We also note that the upper 
critical field H C 2 increases with the increasing 7. 

We now studied the self-consistent results of the order 
parameter and magnetization with 7 = 1.05 in Fig. 2. 
The order parameters are plotted in Figs.2(a)-2(e). As 
seen, for the weaker magnetic field, the order parameter 
is uniform [Fig. 2(a)]. The SC order forms the stripe 
pattern as the exchange field increases to 0.2, as seen in 
Fig. 2(b). The periodicity of the order parameter is about 
24 along y direction as h = 0.20, which is consistent 
with the previous results^ For the case of h — 0.22 and 
h = 0.225, a crossover region from stripe I state to stripe 
II state occurs, indicated in Figs. 2(c) and 2(d). As 
seen in Fig. 2(c), the gap structure has some 2D features 
around the point 'A' with the stripe II state showing 
up gradually. For h = 0.225 [Fig. 2(d)], the stripe II 
state dominates over the stripe I one and the 2D FFLO 
features are almost suppressed. As h increases further, 
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FIG. 2: (Color online) Plots of the order parameter A and 
and the magnetization rm as a function of position for var- 
ious Zeeman fields h with T = 10~ 5 at anisotropic hopping 
parameter 7 = 1.05. 



i.e. for h — 0.23, the SC order forms the stripe pattern 
again with the stripe parallelling to the y axis. 

The spatial distributions of the magnetization (with 
the units %) are shown in Figs. 2(f)-2(j). As seen in 
Fig. 2(f), in the uniform phase, the distribution is also 
uniform, where the magnetization is very weak (about 
0.0017) due to the suppression by the SC order. In the 
ID FFLO stripe I state, as seen in Fig. 2(g), the pat- 
tern also forms ID stripe but the periodicity is only half 
of that of the order parameter. The intensity is largest 
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FIG. 3: (Color online) The LDOS spectra for ID stripe I 
state at h = 0.17. Panel (a) and (b) are the LDOS spectra 
at the nodal line and at the site where the order parameter is 
maximum, respectively. The (blue) dashed and (red) dotted 
lines, and the (black) solid line are spin-up LDOS, spin-down 
LDOS, and whole LDOS, respectively. 



along the nodal lines and is suppressed when the SC or- 
der parameter increases. It reaches the minimum value 
as the SC order is maximum. These features are similar 
to previous results in tetragonal lattice system^. In the 
crossover region, such ash = 0.22, the magnetization also 
behaves the coexistence of the ID and 2D character. As 
h increasing further, the magnetization forms ID stripe 
again with the orientation changing from y axis to x axis. 

Now let us study the LDOS spectra. The LDOS from 
Eq.6 includes two terms, i.e., the spin-up LDOS and spin- 
down LDOS, respectively. In presence of the Zeeman 
field, due to the Pauli paramagnetic effect, the spin-up 
LDOS and spin-down one sperate with the spin-up one 
shifting to the left and the spin-down one to the right, re- 
spectively. The calculated LDOS spectra in the uniform 
phase are similar to the previous report^ which is not 
shown here. The LDOS spectra in the ID stripe I phase 
with h — 0.17 are shown in Fig. 3. Fig. 3(a) is for the 
site on the nodal line. As seen, the spin-up LDOS spec- 
tra show two low-energy peaks locating at u = —0.17 and 
—0.32. The SC coherent peak is suppressed. The spin- 
down LDOS shifts to the right with the two low-energy 
peaks at w = 0.02 and 0.17. These in-gap peaks origi- 
nate from the finite energy Andreev bound states due to 
the sign change in the order parameter across the nodal 
lines, which is the signature of ID FFLO state, as also 
discussed previousl y 22 ' 26 . The peaks at the negative en- 
ergy in total LDOS comes from the spin-up LDOS and 
the peaks at the positive energy are contributed by the 
spin-down LDOS. The intensity of the in-gap peaks will 



FIG. 4: (Color online) The LDOS spectra in the crossover 
region with h — 0.22. Panel (a) is the spectrum at the saddle 
ponit where two nodal lines intersect. Panel (b) is the spec- 
trum at the site where the order parameter is maximum in the 
2D FFLO region. Panels (c) and (d) are the spectra at the 
nodal line and at the site where the order parameter is maxi- 
mum in the ID FFLO region, respectively. The (blue) dashed 
and (red) dotted lines, and the (black) solid line are spin-up 
LDOS, spin-down LDOS, and whole LDOS, respectively. 



decrease as the site moves away from the nodal line. As 
seen from Fig. 3(b), at the site where the order param- 
eter is maximum, the in-gap peaks are turned to be a 
hump at the midgap position for both spin-up and spin- 
down LDOS spectra. The SC coherent peaks are seen 
clearly. The midgap hump is so weak that it is concealed 
in the whole LDOS spectrum. 

The LDOS spectra for the crossover region with h = 
0.22 are presented in Fig. 4. The spectrum at the sad- 
dle point in the local 2D FFLO region [near the point 
"A" in Fig. 2(c)] is presented in Fig. 4(a). Although 
the gap structure is quite different from previous results 
in tetragonal lattice system. While actually for this case 
the spectrum is similar to that of the saddle point's spec- 
trum in the 2D checkerboard FFLO stat o 22 i 26 , namely, 
two kinds of Andreev bound states exist. As a result, 
four in-gap peaks exist in the spin-up LDOS spectra at 
the energies -0.300, -0.208, -0.052, and 0.081 (indicated 
by the arrows). At the site where the order parameter is 
maximum, the LDOS spectrum [Fig. 4(b)] is very com- 
plicated, while the two coherence peaks outside can be 
seen clearly. Shown in Fig. 4(c) is the LDOS spectrum 
on the nodal line in the ID FFLO region. There are 
also two in-gap peaks in both the spin-up and spin-down 
LDOS spectra, indicating the ID characteristics. The SC 
coherent peaks are suppressed again, similar to the case 



FIG. 5: (Color online) Similar to Fig. 3, but for h=0.24. 



of ID stripe phase [Fig. 3(a)]. At the site where the or- 
der parameter is maximum [Fig. 4(d)], the in-gap peaks 
are turned to be a hump at the midgap position for both 
the spin-up and spin-down LDOS spectra, while the SC 
coherence peaks are clearly seen, which leads to the four 
peaks in the whole LDOS spectrum. Summing up the 
above results, it is found that in the crossover region, the 
FFLO indeed includes both the ID and 2D FFLO char- 
acters, which is consistent with the spatial distribution 
of the order parameter and the magnetization (Fig. 2). 

The LDOS spectra of the ID stripe II state with h = 
0.24 is presented in Fig. 5. The result at the nodal line 
is shown in Fig. 5(a). It can be seen that there are two 
in-gap peaks of both the spin-up and spin-down LDOS, 
while the SC coherence peaks are suppressed. Presented 
in Fig. 5(b) is the LDOS spectrum at the site where 
the order parameter is maximum. It is found that the 
in-gap peaks are almost suppressed completely, while the 
SC coherence peaks are clearly seen. Thus, the feature 
of the LDOS spectra of the ID stripe II state is similar 
to that of the ID stripe I state. 

We have shown the LDOS spectra of the three differ- 
ent phases, namely the ID strip 1 phase, the crossover 
region, and the ID stripe II phase. As seen in Fig. 4(a)- 
4(d), the spectra of the Q2D state are quite different 
from that of the other two ID stripe phases. In the Q2D 
state, the ID and 2D FFLO states coexist in different 
regions. These features are expected to be detected by 
the scanning tunneling microscopy (STM) experiments. 

To distinguish the two stripe phases, we present the 
LDOS images with lj = 0. Shown in Fig. 6(a) is the 
LDOS image at h = 0.20 in the ID stripe I state. As seen, 
in this state it also forms ID stripe-like pattern similar to 
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FIG. 6: (Color online) The LDOS maps for quasiparticles 
with w = and different exchange fields. 



the spatial distribution of the order parameter but the pe- 
riodicity is only half of that of the order parameter. The 
intensity is largest along the nodal lines and is suppressed 
when the SC order parameter increases, ft reaches the 
minimum value as the SC order is maximum. The stripe 
parallels to the x axis, consistent with the gap structure. 
The images in the crossover region are presented in Figs. 
6(b) and 6(c) with h = 0.22 and h = 0.225, respectively. 
As seen from Fig. 6(b), the LDOS image in this state 



also behaves the Q2D characteristics, i.e., the stripe- like 
and the local 2D patterns coexist. The intensity is largest 
at the saddle points where two nodal lines intersect. It 
reaches the minimum value as the SC order is maximum. 
The LDOS image at h = 0.23 in the ID stripe II phase is 
presented in Fig. 6(d). It is found that the image forms 
stripe-like pattern again with the orientation switching 
to the y axis. Therefore, it can be used to distinguish 
the two stripe phases by the STM experiments. 

In summary, based on the BdG equations and the 
d-wave superconductivity, the phase diagram and the 
order-parameter structure is studied in presence of the 
external exchange field of a 2D orthorhombic lattice su- 
perconductor. It is found that for weak anisotropic hop- 
ing, i.e. 7 = 1.05, the 2D FFLO state is suppressed and 
only the ID state exists. The orientation of the stripe 
changes as the exchange field increases. The local Q2D 
state exists in the crossover region. These results are 
much different from those of the tetragonal lattice sys- 
tem. Our result suggests that the symmetry of the lat- 
tice plays an important role in the gap structure of the 
FFLO state. The local Q2D state can be detected by the 
STM experiments according to the calculated spectra of 
LDOS. The LDOS images with ui = in these phases 
are also presented, which can be used to distinguish the 
different phases. 
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